Bouncing localized structures in a Liquid-Crystal-Light- Valve experiment 
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Experimental evidence of bouncing localized structures in a nonlinear optical system is reported. 
Oscillations in the position of the localized states are described by a consistent amplitude equation, 
which we call the Lifshitz normal form equation, in analogy with phase transitions. Localized 
structures are shown to arise close to the Lifshtiz point, where non- variational terms drive the 
dynamics into complex and oscillatory behaviors. 
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During the last years localized structures have been 
observed in different fields, such as domains in magnetic 
materials 1], chiral bubbles in liquid crystals 2], cur- 
rent filaments in gas discharge experiments 0, spots in 
chemical reactions 4], pulses 5], kinks |£| and localized 
2D states pj in fluid surface waves, oscillons in granular 
media 8] , isolated states in thermal convection |£J, , 
solitary waves in nonlinear optics 0, 0, El an d cavity 
solitons in lasers [l4| . Localized states are patterns which 
extend only over a small portion of a spatially extended 
and homogeneous system [T^j. Different mechanisms 
leading to stable localization have been proposed 
Among these, two main classes of localized structures 
have to be distinguished, namely those localized struc- 
tures arising as solutions of a quintic Swift-Hohenberg 
like equation 16] and those that are stabilized by nonva- 
riational terms in the subcritical Ginzburg-Landau equa- 
tion 17]. The main difference between the two cases 
is that the first-type localized structures have a char- 
acteristic size which is fixed by the pinning mechanism 
over the underlying pattern or by spatial damped oscil- 
lations between homogenous states 16, whereas the 
second-type ones have no intrinsic spatial length, their 
size being selected by non-variational effects and going 
to infinity when dissipation goes to zero. In both cases, 
non-variational effects may lead to dynamical behaviors 
of localized structures 19]. Variational models based on 
a generalized Swift-Hohenberg equation have been pro- 
posed to describe the appearance of localized structures 
in nonlinear optics [2^. However, a generalization in- 
cluding non-variational terms is generically expected to 
apply even in optics, as ha ppe ns, for instance, in semi- 
conductor laser instabilities [2l[ , giving rise to dynamical 
behaviors of localized structures, such as propagation and 
oscillations of their positions |22j . 

We report here an experimental evidence of local- 
ized structures dynamics in a Liquid- Crystal-Light- Valve 
(LCLV) with optical feedback. It is already known that, 
in the simultaneous presence of bistability and pattern 
forming diffractive feedback, the LCLV system shows lo- 
calized structures 0, HE Q • Recently, rotation of 
localized structures along concentric rings have been re- 



ported in the case of a rotation angle introduced in the 
feedback loop [2£|. Here, we fix a zero rotation angle 
and we show a new dynamical behavior, the bouncing 
of two adjacent localized structures, that is not related 
to imposed boundary conditions but is instead a direct 
consequence of the non-variational character of the sys- 
tem under study. Theoretically, we show that the LCLV 
system has several branches of bistability connecting an 
homogeneous state to a patterned one and we derive an 
amplitude equation accounting for the appearance of lo- 
calized structures. This is a one-dimensional model, that 
we call the Lifshitz normal form equation |22| , character- 
izing the dynamics of localized structures close to each 
point of nascent bistability. 




FIG. 1: Experimental setup: n liquid crystal nematic director; 
Pin and Pfb input and feedback polarizers; L\ and L2 confocal 
25 cm focal length lenses. — L is the free propagation length, 
negative with respect to the plane on which a 1:1 image of 
the front side of the LCLV is formed. 

Description of the experiment. The experimental setup 
is shown in Fig^ The LCLV is composed of a nematic 
liquid crystal film sandwiched in between a glass and a 
photoconductive plate over which a dielectric mirror is 
deposed. The liquid crystal film is planar aligned (ne- 
matic director n parallel to the walls), with a thickness 
d = 15 jam. The liquid crystal filling our LCLV is the 
nematic LC-654, produced by NIOPIK (Moscow) ^ It 
is a mixture of cyano-biphenyls, with a positive dielectric 
anisotropy As = e\\ — e± — +10.7 and large optical bire- 
fringence, An = nil — n± = 0.2, where e\\ and e± are the 
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dielectric permittivities || and _L to n, respectively, and 
n\\ and are n± are the extraordinary (|| to ft) and ordi- 
nary (_L to ft) refractive index, respectively. Transparent 
electrodes over the glass plates permit the application of 
an electrical voltage across the liquid crystal layer. The 
photoconductor behaves like a variable resistance, which 
decreases for increasing illumination. The feedback is ob- 
tained by sending back onto the photoconductor the light 
which has passed through the liquid-crystal layer and has 
been reflected by the dielectric mirror. This light beam 
experiences a phase shift which depends on the liquid 
crystal reorientation and, on its turn, modulates the ef- 
fective voltage that locally applies to the liquid crystals. 

The feedback loop is closed by an optical fiber bundle 
and is designed in such a way that diffraction and polar- 
ization interference are simultaneously present |l2|. The 
optical free propagation length is fixed to L = —10 cm. 
At the linear stage for the pattern formation, a negative 
propagation distance selects the first unstable branch of 
the marginal stability curve, as for a focusing medium. 
The angles of the polarizers are at 45° with respect to the 
liquid crystal director ft. The free end of the fiber bundle 
is mounted on a precision rotation and translation stage, 
to avoid rotation or translation in the feedback loop. 
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FIG. 2: Spatially averaged feedback intensity < I w > (units 
are the gray values, g.v., on the CCD camera) as a function 
of Vb; input intensity I in = 0.75 mW/cm 2 . 

For this parameter setting, as Vb increases there is a 
series of successive branches of bistability between a pe- 
riodic pattern and a homogeneous solution. In FigOwe 
report the spatially averaged feedback intensity < I w > 
measured for a fixed value of the input intensity, Zj n = 
0.75 mW/cm 2 , and for varying Vb, by integrating the 
images on the near-field CCD camera (see Fig[Q . The 
abrupt changes of < I w > correspond to the appearance 
of localized structures and thus roughly indicate the loca- 
tions of the nascent bistability points. The peak value in- 
tensity of the localized structures is approximately twice 
the average value < I w >. In the LCLV system, the 
bistability between homogenous states results from the 
subcritical character of the Freedericksz transition, when 
the local electric field, which applies to the liquid crystals, 
depends on the liquid crystal reorientation angle [23, |2£j . 
Here, we limit our study to the bistable branch located 
around Vb = 13.2 V rrns (frequency 5 KHz), however sim- 
ilar observations can be obtained close to any other of the 



nascent bistability points. 

We have carried out one-dimensional experiments, in 
order to avoid the influence of any optical misalignment 
(such as small drifts) on the dynamics of localized states. 
A rectangular mask is introduced in the optical feedback 
loop, just in contact to the entrance side of the fiber bun- 
dle. The width of the aperture is D = 0.50 mm whereas 
its length is / = 20 mm. The size of each localized struc- 
tures is A ^ 350 /im, so that the transverse aspect ratio 
D/A ~ 1 is small enough for the system to be consid- 
ered as one-dimensional and the longitudinal aspect ratio 
1 1 A ~ 60 is large enough for the system to be considered 
as a spatially extended one. 
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FIG. 3: Snapshots showing bouncing localized structures; a) 
t = 0.0, b) 1.0, c) 1.3, d) 1.7, e) 2.1, f) 2.4, g) 2.8 sec. 




FIG. 4: Localized structures profile, I w {xo,y); xo is the loca- 
tion of the dashed line in Fig^] a) t = 0.0 b) 1.3, c) 1.7 sec. 

Instantaneous snapshots of bouncing localized struc- 
tures, together with their spatial profiles, are shown in 
Figdl for Vb 13.2 V rms and I in = 0.95 mW/cm 2 . 
The corresponding spatio-temporal plot is displayed in 
Figinb, showing the periodic oscillations for the positions 
of the structures. In the same figure, FigJSt, it is shown 
the spatio-temporal plot corresponding to stationary lo- 
calized structures, as observed for a slightly decreased 
input intensity, Ii n = 0.90 mW/cm 2 , and for the same 
value of Vb- Figl^t displays the spatio-temporal diagram 
corresponding to aperiodic oscillations in the structure 
positions, as observed for Vb = 13.3 V rms and Ii n = 0.90 
mW/cm 2 . The dynamical behavior of localized struc- 
tures is very sensitive to parameter changes and, even 
though their appearance is clearly located around each 
point of nascent bistability, their stability range is smaller 
than the width of the bistable region. When loosing sta- 
bility, localized structures either form clusters or anni- 
hilate, depending if they are driven on the pinning or 
depinning side of the bistable region 18]. However, a 
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careful experimental characterization of the pinning front 
is a work still in progress. 




FIG. 5: Space (vertical) -time (horizontal) diagrams showing 
a) two stationary localized structures, b) periodic and c) ape- 
riodic oscillations of the structure positions. The total elapsed 
time is 94 s. 

Theoretical description. The light intensity I w 
reaching the photoconductor is given by I w — 

h^\ e ~ l ik d ^(\ -f e -i(3cos 6>j|2^ w h ere x i s the transverse 
direction of the liquid crystal layer, f3cos 2 is the overall 
phase shift experienced by the light travelling forth and 
back through the liquid crystal layer; (3 = 2kdAn, where 
k = 2tt/X is the optical wave number (A = 633 nm). 

As long as Ii n is sufficiently small, that is, of the or- 
der of a few mW/cm 2 , the effective electric field, E e ff, 
applied to the liquid crystal layer can be expressed as 
E e ff = TVo/d + al w , where Vq is the voltage applied to 
the LCLV, < T < 1 is a transfer factor that depends 
on the electrical impedances of the photoconductor, di- 
electric mirror and liquid crystals and a is a phenomeno- 
logical dimensional parameter that describes the linear 
response of the photoconductor [2^ . 

Let us call 6 (x, t) the average director tilt. 6 = 
is the initial planar alignment whereas = ir/2 is the 
homeotropic alignment corresponding to saturation of 
the molecular reorientation. The dynamics of is de- 
scribed by a local relaxation equation of the form 

with V eff = E eff d = rVb + adl w (0,d x ) > V FT the ef- 
fective voltage applied to the liquid crystals, Vft the 
threshold for the Freedericksz transition and I the elec- 
tric coherence length. The above model have been de- 
duced by fitting the experimental data for the open loop 
response of the LCLV |22| and it is slightly different with 
respect to the one proposed in Ref. |l2|. It is important 
to note that a rigorous derivation of the response func- 
tion of the LCLV would require a modal expansion along 
the longitudinal direction of the liquid crystal layer. 
The homogeneous equilibrium solutions are 6q = 

when V e ff < Vft and Oq = tt/2 (l — ^JVft ~/Vef 7) 
when V e ff > Vft- Above Freedericksz transition and 
by neglecting the spatial terms we can find a closed 
expression for the homogeneous equilibrium solutions: 

#0 = tt/2 (l - y/V FT /(rV + al in [l + cos(/3cos 2 O )])). 



The value of Vft is set to 3.2 V rms , as measured for 
the LCLV HI Q, and the graph of o (V o J in ) is plot- 
ted in FigEl In agreement with the bistability branches 
observed experimentally, several points of nascent bista- 
bility can be distinguished, corresponding to the critical 
points where 6q(Vq, Ii n ) becomes a multi- valued function. 




FIG. 6: The multi- valued function Oo(Vo,Ii n )- Shaded areas 
show the location of the nascent bistability points. 

Close to each point of nascent bistability, and neglect- 
ing spatial derivatives, we can develop 6 = 6q + u + ... and 
derive a normal form equation describing an imperfect 
pitchfork bifurcation [l5[, dtu = 77 + /^ — u 3 + h.o.t, where 
ji is the bifurcation parameter and 77 accounts for the 
asymmetry between the two homogeneous states. Higher 
order terms are ruled out by the scaling analysis, since 
u ~ /i 1 / 2 , rj ~ /i 3//2 and dt ~ /i, [i <C 1. If we now con- 
sider the spatial effects, due to the elasticity of the liquid 
crystal and to the light diffraction, the system exhibits a 
spatial instability as a function of the diffraction length 
and, since the spatial dependence of I w is nonlocal, the 
dynamics is a non- variational one. 

The confluence of bistability and spatial bifurcation 
give rise to a critical point of codimension three, that we 
call the Lifshitz point, in analogy with the triple point in- 
troduced for phase transitions in helicoidal ferromagnetic 
states Close to this point, we derive an amplitude 
equation, that we call the Lifshitz normal form |22j . 

d t u = rj -\- i±u — -\- v d xx u — d xxxx u 

+dud xx u + c(d x ii) 2 , (2) 

where d x ~ /i 1 / 4 , v ~ /i 1 / 2 accounts for the intrinsic 
length of the system (diffusion), d ~ O (1) and c ~ O (1). 
The term d xxxx u describes a super-diffusion, accounting 
for the short distance repulsive interaction, whereas the 
terms proportional to d and c are, respectively, the non- 
linear diffusion and convection. The full and lengthy ex- 
pressions of these coefficients, as a function of the LCLV 
parameters, will be reported elsewhere |2J|. Note that 



4 



the same model has been recently deduced for instabili- 
ties in a semiconductor laser |2l| . 

The model shows bistability between a homogeneous 
and a spatially periodic solutions and therefore exhibits 
a family of localized structures. Depending on the choice 
of parameters, localized structures may show periodic or 
aperiodic oscillations of their position. We can fix \i and 
r] by varying Vq and Ij n . More interesting is the be- 
havior of the effective diffusion term which has the 
form v oc I 2 + (tt(3Lcos 2 ((3/2cos 2 o )sm20 )/4k((TV o + 
aIi n (l+cos(/3 cos 2 Oo))). Only when the optical free prop- 
agation length L is negative it is possible, by increasing 
the input intensity Ij n , to drive the system through the 
Lifshitz point (y changes its sign from positive to neg- 
ative). This means that stable localized structures can 
be obtained only for a focusing Kerr-like nonlinearity. 
Once crossed the Lifshitz point, the nonlinear diffusion 
coefficient is negative and the convection coefficient is 
positive. In this region of parameters, numerical simu- 
lations of Eq.(|2j show a qualitative agreement with the 
experimental observations, as shown in Fig0 
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FIG. 7: Numerical simulations of Eq. J5J for n — —0.02, 
/a -0.02, v -1.00, c 2.00 and d -1.51, showing two 
bouncing localized structures (spatial profile in the inset). 

Conclusions. We have shown a new kind of localized 
structure dynamics, consisting of a bouncing behavior 
between two adjacent structures, and we have described 
it by an universal model, the Lifshitz normal form equa- 
tion. The Lifshitz equation, that reduces to the gen- 
eralized Swift-Hohenberg equation for n = d = c = 0, 
has already been used to describe the transition from 
smectic to helicoidal phase in liquid crystals l32l and the 
pulse dynamics in reaction diffusion systems [33| . When 
one neglects the cubic and the nonlinear diffusion terms 
(d = 0), it reduces to the Nikolaevskii equation that de- 
scribes longitudinal seismic waves jiij . 
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